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ABSTRACT 

Since 1950, when Oort published his hypothesis several important new facts were established 
in the field. In this situation the apparent source region (or regions) of long-period comets as 
well as the definition of the dynamically new comet are still open questions, as well as the 
characteristics of the hypothetical Oort Cloud. The aim of this investigation is to look for the 
apparent source of selected long period comets and to refine the definition of dynamically new 
comets. 

Basing on pure gravitational original orbits all comets studied in this paper were widely called 
dynamically new. We show, however, that incorporation of the non-gravitational forces into 
the orbit determination process significantly changes the situation. 

We determined precise non-gravitational orbits of all investigated comets and next followed 
numerical ly the ir past and future motion during one orbital period. Applying ingenious 
Sitarski's (119981) method of creating swarms of virtual comets compatible with observations, 
we were able to derive the uncertainties of original and future orbital elements, as well as the 
uncertainties of the previous and next perihelion distances. 

We concluded that the past and future evolution of cometary orbits under the Galactic tide 
perturbations is the only way to find which comets are really dynamically new. In our sam- 
ple less than 30% of comets are in fact dynamically new. Most of them had small previous 
perihelion distance. On the other hand, 60% of them will be lost on hyperbolic orbits in the 
future. This evidence suggests that the apparent source of long-period comets investigation is 
a demanding question. 

We also have shown that a significant percentage of long-period comets can visit the zone of 
visibility during at least two or three consecutive perihelion passages. 
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1 INTRODUCTION 



Jan Oortl ll95Ql) proposed a hypothesis, that the Solar System is 
surrounded by a huge, spherical cloud of comets, now called the 
Oort Cloud. As the argument for the existence of such a cloud 
he completed a list of 19 original orbit|3 of well observed long- 
period (hereafter LP) comets and showed, that their inverse of semi- 
major axis 1/fl has the distribution apparently peeked near zero, 
at the positive side. He concluded that new long-period comets 
come from distances of 50000 to 150000 AU. He also showed, 
that perturbations by passing stars can change cometary orbit sig- 
nificantly, making it observable as "dynamically new" long-period 
comet. During the last 60 years several new factors occurred in this 
field. 
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' In this paper we define original and future orbits as barycentric osculating 
orbital elements derived at a distance of 250 AU from the Sun. 



First, the population of precise original cometary orbits in- 
creased to several hundreds. In the rec ently published, 17th edition 
of the Catalogue of Cometary Orbits jMarsden & William"sll2008l) 
(hereafter MWC08) 499 precise original and future cometary orbits 
are listed. Second, new and dominating perturbing force must be 
included in the model: tidal action of the Galactic disk and Galac- 
tic centre. They are negligible when calculating original and future 
orbits but when following the past or future motion of a comet at 
larger distances one must take them into account. It means, that 
looking at the original orbit elements we cannot tell the past or 
predict the future motion of this body, for example foretell the 
previous/ next perihelion or aphelio n distances. This was demon- 
strated bv lDvbczvhskil j200ll : l2006h . Third, we learned how to in- 
vestigate and in many cases we can successfully determine non- 
gravitational forces in the motion of long-period comets. This can 
significantly chan ge our knowledge of original and future orbit s 
(see for example: iMarsden et alll 19731 : lKr61ikowskall200lL l2006h . 
It also means, that tens of precise original orbits with determined 
non-gravitational forces could be added to the MWC08 list. 



2 M. Krdlikowska and P.A. Dybczynski 



Taking all this into account, the apparent source region (or 
regions) of long-period comets as well as the definition of the dy- 
namically new comet are still open questions, as well as the char- 
acteristics (if not the existence) of the hypothetical Oort Cloud. 
In this paper we start to investigate these problems concentrat- 
ing on 26 comets constituting so-called non-gravitational (hereafter 
NG) Oort spike, i.e. having l/aojiNc < lO^'* AU^'. This particu- 
lar choice comes from the fact, that after applying non-gravitational 
model in orbit determination we obtain significantly different dis- 
tribution of the inverse original semi-major axes. Almost all 1 /oori 
values were shifted to larger values and among them a lot of hyper- 
bolic cases changed into elliptic ones. Since we want to study the 
implications of obtaining non-gravitational original orbits on the 
Oort hypothesis, the past (and future) motion of these 26 comets 
were analyzed at first. 

In Section |2] we describe observational material and its pro- 
cessing. In the next section calculation of original and future 
cometary orbits for the selected 26 non-gravitational Oort spike 
comets is described in detail. Section|4]describes the numerical in- 
tegration of the past and future motion of these comets under the ac- 
tion of Galactic tides. Then, results and conclusions are presented. 



2 OBSERVATIONS AND THEIR PROCESSING 

MWC08 includes 184 comets with 1 /nori < 10"'* AU"' where for 
majority of them the orbits are determined with the highest ac- 
curacy. At first we looked at 154 most precise orbits. Thus, we 
took into account all comets with the qu ality class 1 according to 
the classification introduced by Marsden et al.l ( Il978l) . About half 
of them have perihelion distances less than 3.0 AU, what means 
that we should suspect detectable influence of the non-gravitational 
forces. The inverse dependence of the strength of non-gravitational 
forces on the perihelion distance is widely expected and indeed is 
clearly visible in our Fig|5](to be described in detail in the next 
section). We took into account almost all these objects (except 
ten comets discovered before 1950 where positional data were not 
available for us) and additionally all LP comets with the NG orbits 
in MWC08 (29 objects, 1 /aoi-[ are not given for them in the Cat- 
alogue). In majority of cases, the observational material taken for 
orbital determinations is more complete than in the MWC08: for 
comets from the XIX century the observations were collected di- 
rectly from the published papers, for comets from the first half of 
XX century the data were collected in Warsaw in cooperation with 
the Slovakian group at the Astronomical Institute in Bratislava and 
Tatranska Lomnica, more modern data were taken from the Mi- 
nor Planet Electronic Circulars available through the Web Page at 
|http: //www.cfa.harvard.edu 

^From these class 1 objects we were able to derive the sam- 
ple of 26 LP comets with determinable non-gravitational effects 
from the astrometric data and with the non-gravitational recipro- 
cals of semi-major axes not greater than lO"'* AU"' (l/flori,NG < 
10"^ AU"' , hereafter: NG Oort spike). Among them only 9 comets 
(C/1885 XI, C/1956 Rl, C/1959 Yl, C/1986 PIA, C/1989 Ql, 
C/1990 Kl, C/1993 Al, C/2001 Q4, C/2002 T7) have NG orbits 
presented in MWC08. C/1990 Kl Levy has formally I/ooiing 
slightly greater than the adopted limit but the uncertainty inter- 
val of it evidently overlap with the NG Oort spike defined above. 
We have not included in our sample objects which are still poten- 
tially observable because the derived orbits for these objects could 
not be treated as definitive results. The preliminary investigation 
of the original non-gravitational reciprocals of semi-major axes. 



Table 1. Observational material for 26 comets under consideration 



No. 

Name Observational arc obs. 



1 /oori,NG> given in lKrolikowskaldioOdh for the majority of these 
26 comets. In that paper data have not been weighted, while now, to 
obtain the more reliable orbits, we decided to include the weighting 
of observations into the data processing. 

The determination of the NG parameters in the motion of 
LP comets is very difficult mainly due to limited observational 
material covering one apparition or even just half apparition (in 
the case when comet was discovered after its perihelion passage). 
Thus, the processing of astrometric data is crucial for this purpose. 
iKiolikowsk a ( 2006|) divided each set of astrometric data into pre- 
perihelion and post-perihelion arc (when it was possible) and then 
selected data in each part independently. The difference of rms 
for the whole data set (for the pure gravitational solution) and the 
rms obtained separately for both parts of the orbit indicates how 
strongly the NG-effects affect the orbital motion. After many tests 
we have recognized, however, that the weighting is crucial for the 
orbit fitting for comets discovered before 1950 and sometimes for 
modern comets, too. Thus, we decided to adopt more advanced data 
treatment. Each set of astrometric data has been processed indi- 
vidually for the pure gravitational orbit and non-gravitational case 
independently in the following way. 

For each comet we tested two methods of selecting and 
weighting data for both: pure gravitational and non-gravitational 
orbit, independently. It means that we determined four nominal so- 
lutions for each object: two for pure gravitational case and two for 
non-gravitational case (based on different data treatment). It allows 
us to choose the best and as homogeneous as possible data process- 
ing for the whole sample of LP comets. We tried to use the most 
but reasonable discriminative selection criterion and according to 
our tests the year 1950 can serve as a suitable dividing date. Thus, 
we divided all considered comets into two sets: comets discovered 
before and after 1950. These two sets of observations were refined 



C/1885 XI Fabry 
C/1892Q1 Brooks 
C/1913 Yl Delavan 
C/1940 R2 Cunningham 
C/1946 Ul Bester 
C/1952 Wl Mrkos 
C/1956 Rl Aiend-Roland 
C/1959 YlBurnham 
C/1978 HI Meier 
C/1986 PI Wilson 
C/1989 Ql Okazaki-Levy-R. 
C/1989 XI Austin 
C/1990 Kl Levy 
C/1991 F2 Helin-Lawrence 
C/1993 Al Mueller 
C/1993 Ql Mueller 
C/1996 El NEAT 
C/1997 BA6 Spacewatch 
C/1997 J2 Meunier-Dupouy 
C/1999Y1 LINEAR 
C/2000 SV74 LINEAR 
C/2001 Q4NEAT 
C/2002 E2 Snyder-Murakami 
C/2002 T7 LINEAR 
C/2003 K4 LINEAR 
C/2004B1 LINEAR 



1885 12 02 - 


1886 07 20 


228 


1892 09 01 - 


1893 07 13 


191 


1913 10 25 - 


1915 09 07 


1009 


1940 08 25 - 


1941 04 01 


370 


1946 11 01 - 


1948 10 02 


142 


1952 12 10- 


1953 07 18 


36 


1956 11 08 - 


1958 04 11 


249 


1960 01 04- 


1960 06 17 


88 


1978 04 28 - 


1979 12 09 


287 


1986 08 05 - 


1989 04 11 


688 


1989 08 24 - 


1989 12 24 


231 


1989 12 06- 


1990 06 27 


281 


1990 05 21 - 


1992 04 01 


678 


1991 02 23- 


1992 09 30 


114 


1992 11 26- 


1994 08 17 


746 


1993 08 16- 


1994 04 17 


526 


1996 03 15 - 


1996 10 12 


249 


1997 01 11 - 


2004 09 15 


529 


1997 05 05- 


1999 10 09 


1446 


1999 10 29- 


2003 07 19 


884 


2000 09 05 - 


2005 05 12 


2189 


2001 08 24- 


2006 08 18 


2661 


2002 03 08 - 


2003 01 08 


941 


2002 10 12- 


2006 03 07 


4451 


2003 05 28 - 


2006 11 17 


3658 


2004 01 28 - 


2008 08 24 


2057 
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Figure 1. The O-C distributions for Comet C/2004 Bl LINEAR and C/2()01 Q4 NEAT. Tlie non-weiglited data are presented in the left-side column for 
each comet and the weighted data - in the right column where the gravitational cases are displayed in upper figures and the NG cases in lower figures. The 
best-fitting Gaussian distributions are shown by dots, each chosen exactly in the middle of histogram bin. 



independently for gravitational and non-gravitational orbit determi- 
nation by applying 

for comets before 1 950: 

• Chauvenet's criterior0 for selection procedure + weighting, 

• Bessel's criterion^ for selection procedure + weighting, 

for comets after 1 950: 

• independent Bessel's selection for observations before and af- 
ter perihelion, no weighting, 

• Bessel's criterion for selection procedure + weighting. 

To reduce systematic errors in the observational material, such as 
the bias associated with a site as a function of time, one should 
consider specific procedures. For example, in the case of long se- 
quences of observations we divided the whole observational mate- 
rial into a few time subintervals according to the internal structure 
of material, i.e. according to the existing gaps in observations. 

The data of two comets with the number of observations less 
than 100 (C/1952 Wl and C/1959 Yl) have not been weighted. 
For two of the comets discovered before 1950 the Chauvenet's and 
Bessel's criterions gave the same set of residuals and for the next 
three comets the differences in the number of rejected residuals are 
very small (less than 1%). Thus, we present here results derived for 
more discriminative Bessel's criterion also for comets before 1950. 
One should note that our data processing resulted in reasonably 
low number of rejections (below 17% - for comets before 1960 
and with mean percentage of rejected data of 3.4% - for comets 
after 1960; the only extreme case here is C/2004 Bl, for which 9% 
of observations was rejected). The description of observational data 
used for orbital determination is presented in Table[T]where comets 
are ordered by discovery date. 



short description of both criterions is given in l7.1l 



To determine th e NG cometary orbit w e used the standard for- 
malism proposed bv lMarsden etal j ( ll973h where the three orbital 
components of the NG acceleration acting on a comet are symmet- 
ric relative to perihelion: 



Fi=Ai- g{r) 



A[ = const for i = 
-a{r/ro)-"'[l + ir/roy 



(1) 

(2) 



where Fj , F2, F3 represent the radial, transverse and normal com- 
ponents of the NG acceleration, respectively. The exponential co- 
efficients m, n, k are equal to 2.15, 5.093, and 4.6142, respectively; 
the normalization constant a = 0.1113 gives g(l AU) = 1; the 
scale distance ro = 2.808 AU. From orbital calculations, the NG pa- 
rameters Aj ,A2, and A3 were derived together with six orbital ele- 
ments within a giv en time interval (numerical details are given in 
lKr61ikowsk3l2006l) . This NG model assumes that water sublimates 
from the whole surface of an isothermal cometary nucleus. One can 
see in the 10* column of Table|2]that in some cases just two NG pa- 
rameters were determinable with reasonable accuracy (in two cases 
- only one NG parameter). We decided to use in this investigation 
the minimal necessary number of NG parameters. For example, in 
the case of C/1993 Al and C/1996 El the NG parameter A3 is de- 
terminable, however, the decrease of rms (relative to NG solution 
with A 1 and A2) is insignificant with no perceptible changes in O- 
C diagrams and O-C distributions (compare with online material in 
lKr61ikowsk3 l l2004l) ). The asymmetric NG model which introduces 
the additional NG parameter t - the time displacement of the max- 
imum of the g {r(t — t)) relative to perihelion, seems be important 
just for one comet, C/1959 Yl Burnham, what makes the number 
of NG parameters for this comet to be equal 4. 

To decide whether the derived NG solution better represents 
the actual motion of individual comet than the pure gravitational 
orbit we use three criterions; 
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C/1997 J2 (GR solution: weighted data) 
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• decrease of the rms for the NG solution in comparison with 
the pure gravitational solution, 

• smaller trends in 0-C diagrams for NG orbit fitting, 

• the 0-C distribution of NG orbit closer to the normal one. 



1999.8 After inspection of the O-C diagrams and distributions we always 
conclude that the NG orbit derived from the weighted and selected 
data sets represents the actual comet motion much better than the 
pure gravitational orbit. This fact is visualized in Fig. [T] for two 
comets C/2004 Bl and C/2001 Q4. Therefore, the pure gravita- 
tional orbits are presented in this paper only for comparison and 
to show that such orbits could produce misleading orbital results 
and incorrect conclusions about future and past history of cometary 
dynamics. The examples of differences in the 0-C diagrams be- 
tween pure gravitational and NG solutions are presented for two 
comets in Figs. |2] and [S] For comet C/1997 J2 only pure gravita- 
tional orbit is presented in MWC08 as well as at the JPL Small- 




Figure 2. The O-C diagrams for Comet C/1997 J2 Meunier-Dupouy 
(weighted data). Two upper figures are given for O-C based on pure gravi- 
tational solution, two lower - for 0-C based on NG solution. Residuals in 
right ascension are shown as black dots and in declination - as grey crosses; 
the moment of perihelion passage is shown by dashed vertical line. 
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Figure 3. The O-C diagrams for Comet C/2001 Q4 NEAT (weighted data). 
Two upper figures are given for O-C based on pure gravitational solution, 
two lower - for O-C based on NG solution (for the clarity the three last 
observations taken in August 2006 are not displayed). Symbols have the 
same meaning as in Fig.|2] 



Body Database Browser ( ht tp://ssd.jpl.nasa. gov/sbdb.cgi^ , whereas 
for Comet C/2001 Q4 NG solutions are presented in both sources. 
Fig |2] reveals that some trends visible in pure gravitational 0-C 
diagram for C/1997 J2 Meunier-Dupouy are almost canceled for 
NG solution while the rms decreased from 0"91 to 0!'67. Similar 
behavior was noticed for all NG solutions given in this paper even 
for four objects with largest perihelion distance (q> 3.0 AU) as was 
shown in Fig.|2]for Comet C/1997 J2. 

Among investigated sample of NG Oort spike comets 
C/2001 Q4 NEAT presents the case of one of the most sig- 
nificant decrease of rms from l"69 (pure gravitational orbit, 
non-weighted data) to 0!'63 (NG orbit, weighted data) (Ta- 
ble |2l Figs. Q] [3j. In this group are also C/1959 Yl Bumham, 
C/1990 Kl Levy, C/1993 Al Mueller and C/2002 T7 LINEAR. 

For many comets with evidently detectable NG-effects we ob- 
tained strongly non-Gaussian O-C distribution for pure gravita- 
tional orbit whereas the O-C distribution follows the normal dis- 
tribution when the NG orbit were fitted. An excellent example of 
such behavior is Comet C/2004 Bl. The departures from the best 
Gaussian fitting to 0-C residuals are shown in Fig.[T] One can note 
that the deviations of 0-C distribution from normal distribution in 
the pure gravitational case always decreased for weighted data in 
comparison to non-weighted data, however in some cases this devi- 
ation is still clearly visible even for weighted data as one can see for 
both comets presented in Fig.[T] It seems evident that non-Gaussian 
distribution of 0-C residuals results from ignoring of NG forces 
in the cometary motion. In a few cases (C/1990 Kl, C/1993 Al, 
C/2001 Q4, C/2002 T7 and C/2003 K4), however, the 0-C distribu- 
tions for NG solution obtained on the basis on weighted data differ 
from the Gaussian distribution also. In all these five cases the resid- 
uals in NG model still display some trends too (especially around 
perihelion passage), though less significant than in the pure gravi- 
tational case (Fig.O. We tested that the time-shifted g(r)-function 
relative to perihelion passage do not improved the O-C diagrams 
here. It seems that some other function than standard g(r) would 
be more adequate to describe the NG effects in these comets. It 
should be noticed, however, that the presented NG nominal orbits 
for these five comets are closer to actual cometary orbit than pure 
gravitational solutions. 

And the last, important remark concerning the gravitational 
solutions determined from weighted data. In practice, these so- 
lutions are not purely gravitational since the time distribution of 
weights could compensate the NG trends. For three objects from 
our sample (C/2002 T7, C/2001 Q4 and C/1997 BAe) relative de- 
creasing of rms and some time-dependence in weights are evidently 
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visible . For the most spectacular case (C/2002 T7) the rms of 
gravitational solutions decreased from 1 "21 (non-weighted data) to 
0"61 (weighted data) whereas the rms of NG solutions decreased 
from 0"64 to 0"58, respectively, and systematic trends in weights 
disappeared. Thus, we should keep in mind that so called gravi- 
tational solutions determined on the basis of weighted data could 
partially incorporate NG trends. 
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3 OBTAINING THE MOST PROBABLE ORIGINAL AND 
FUTURE ORBITS 

To investigate the past and future history of each particular comet 
we have examined the evolution of thousands of virtual comets 
(VCs) orbits from the confidence region, i.e. the 6-dimensional re- 
gion (7-9D region in the NG case) of orbital elements (and NG pa- 
rameters), each time fully compatible with the observations. We 
constructed the confidence region using the Sitarski's method of 
the random orbit selection l lSitarskilll998l) . 

This method allows us to generate any number of randomly 
selected orbits of VC. It should be noticed that according to this 
random selection method the derived sample of VCs follows the 
normal distribution in the orbital elements space (also the rms's ful- 
fil the 6-9 dimensional normal statistics; Rrolikowska et al. 200^. 

In the calculations described here we fill a confidence region 
with 5 000 VCs for each nominal solution (GR and NG) for each 
considered comet. Fig.|4]shows projections of the 6D/8D parameter 
space of 5 000 VCs of C/1989 XI (upper panels) and C/1993 Al 
(lower panels) onto the plane of two chosen orbital elements. Or- 
bital cloning procedure was applied close to the observational arc 
(at the epoch of 1990 April 19 and 1992 October 25, respectively). 
The derived swarm of VCs follows the normal distribution in the 
6D space of orbital elements (8D space of orbital elements and two 
NG parameters in the NG case for both comets). This is visual- 
ized by four grey tints of points in Fig. [4] Each point represents 
a single VC, while its grey tint indicates the deviation magnitude 
from the nominal orbit with the confidence levels defined as fol- 
lows: < 50 %, 50 % - 90 %, 90 % - 99 %, and > 99 % (from the 
darkest to brightest points, respectively). The symbols in crowded 
areas heavily overlap and the brighter points are often overprinted 
with darker ones. One can see that the NG swarm of C/1989 XI 
is more disperse than its pure gravitational swarm (upper panel 
of Fig. |4l(. Similar situation took place in almost all investigated 
comets except C/1993 Al and C/2000 SV74. The largest dispersion 
of original NG swarm within our sample of NG Oort spike comets 
appeared in case of comet C/1959 Yl Bumham where four NG pa- 
rameters were determined. In the C/1993 Al case the compactness 
of NG swarm of VCs is higher than compactness of GR swarm 
though in the first case we determined two more parameters from 
the same data set than in the latter case (lower panel of Fig.l?^. It is a 
consequence of radically better NG orbit fitting to data (rms=0"98) 
in comparison to pure gravitational orbit fitting (rms=2"49). In the 
majority of cases, however, the relative decreasing of rms is smaller 
and the NG swarm of VCs is more disperse than pure gravitational 
swarm, as expected when increasing parameter space dimensions. 

To calculate the original and future swarm of cometary or- 
bit each of VC was followed from its position at osculation 
epoch backwards and forwards until the VC reached a distance of 
250 AU from the Sun. The equations of comet's motion have been 
integrated numerical ly using the recurrent power series method 
( ISitarski|[l989l . l2002h . taking into account perturbations by all the 
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Figure 4. Projection of the 6D/8D space of 5 000 VCs of C/1989 Ql (upper 
panel) and C/1993 Al (lower panel) onto the e-q plane. Left plots show 
the pure gravitational swaiTns of VCs, right plots - NG swarms. Each point 
represents a single virtual orbit, while its grey tint indicates the deviation 
magnitude from the nominal orbit with the confidence level of: < 50 %, 
50 % - 90 %, 90 % - 99 %, and > 99 % (from the most dark grey points 
to the most light grey points, respectively). Each plot is centered on the 
nominal values of respective pair of osculating orbital elements denoted by 
the subscript '0'. 



planets and including the relativistic effects. In this way we were 
able to obtain the nominal original/future orbit of each comet as 
well as the uncertainties of the derived values of orbital elements 
by fitti ng the normal distrib ution to each original/future cometary 
swarm l lKr61ikowsk3l200ll) . Original and future semi-major axes 
with their uncertainties are given in Table |2] where comets are or- 
dered by discovery date. Fig. |5] visualizes the differences between 
original 1/a derived for NG and GR swarms of comets. One can 
note that all investigated comets have l/aori,GR < 5 x 10^^ AU^' . 
This means, t hat al l these comets would belong to the first row of 
Table 1 in the lOorll l ll950l) paper. 



4 PAST AND FUTURE MOTION UNDER THE 
GALACTIC TIDE 

4.1 Force model 

FoUowing lFouchard et alj J2005h and others, we used a set of equa- 
tions of motion in the vicinity of the Sun, accounting for Galactic 
disk tide and Galactic centre tide. Pretka ( 1998, 1999) argued that 
for relatively small time intervals (say ten million years) it is not 
necessary to include the Galactic centre term, as its influence man- 
ifests significantly on considerably longer time intervals. However, 
to maintain compatibility with other works in the field and taking 
into account, that including Galactic centre term into the model of 
Galactic tides is not very time consuming we decided to use such 
a full model. During all integrations we performed calculations in 
both: simplified and full Galactic tide models, just to confirm the 
above-mentioned conclusions and observe the difference in some 
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Table 2. Original and future semi-major axes derived from pure gravitational nominal solutions (columns 2-3) and NG nominal solutions (columns 4-5) where 
the number of NG parameters determined for NG solutions is given in the column 10. The rms's and number of residuals are given in the columns 6-7 and 
8-9, respectively, where 'GR' refers to gravitational nominal solution. If there are two rows for a particular comet, the first describes results with weighting 
and the second without weighting. Single row (older comets) means that weighting was applied. 



Name gravitational solutions NG solutions rmsQR rmsfjo number of 





l/aori 




1/afut 




1/aori 




1/afut 












NG 




in units 


of 10"^ AU 


-1 


in units 


of 10"'' AU 


-1 


II 


// 






par. 


1 


2 




3 




4 




5 




6 


7 


8 


9 


10 


C/1885 XI Fabry 


-2.73± 11.28 


-256.73± 11.28 


-l-61.04± 


17.17 


-103.42± 


41.33 


3.92 


3.58 


394 


390 


2 


C/1892 Ql Brooks 


-57.34± 12.52 


-568.73± 12.52 


-l-58.07± 


50.88 


-486. 86± 


23.28 


3.20 


2.80 


343 


334 


2 


C/1913 Yl Delavan 


-i-22.14± 


1.30 


+56.41± 


1.30 


-h52.57± 


4.23 


+86.84± 


4.23 


2.06 


2.00 


1874 


1860 


1 


C/1940 R2 Cunningham 


-60.00± 17.66 


-1630.62± 17.67 


-i-51.66± 


17.30 


-1494.22± 


36.52 


1.78 


1.49 


678 


670 


2 


C/1946 Ul Bester 


-l-5.72± 


4.97 


+32.53± 


4.97 


-h25.41± 


5.25 


-l-52.21± 


5.25 


1.53 


1.35 


257 


253 


1 


C/1952 Wl Mrkos 


-140.69± 21.69 


-298.66± 21.69 


-0.07± 


85.78 


-41.39± 106.94 


1.14 


1.08 


61 


61 


2 


C/1956 Rl Arend-Roland 


-97.97± 


5.46 


-603.5 1± 


5.47 


-Hl0.37± 


11.42 


-579.69± 


9.97 


1.68 


1.39 


461 


458 


3 




-104.38± 


5.02 


-609.92± 


5.03 


-l-7.83± 


11.19 


-587.32± 


10.59 


1.68 


1.41 


449 


449 


3 


C/1959 Yl Burnham 


-138.32± 15.93 


-589.97± 15.94 


-l-2.09± 123.88 


-286.65± 


33.46 


3.00 


1.60 


146 


146 


4 


C/ 1978 HI Meier 


+23.79± 


2.42 


-1028.68± 


2.42 


-h80.55± 


10.73 


-1025. 19± 


8.15 


1.04 


1.00 


565 


565 


2 




+24.52± 


2.73 


-1027.96± 


2.73 


-h89.23± 


12.83 


-1018.74± 


8.09 


1.22 


1.18 


549 


549 


2 


C/1986 PI Wilson 


-i-9.53± 


0.82 


-1-735. 14± 


0.82 


-h42.00± 


1.97 


-h767.22± 


2.35 


1.35 


1.11 


1364 


1361 


2 




-i-0.44± 


0.99 


-i-725.17± 


0.99 


-h42.33± 


2.39 


-h767.22± 


2.86 


1.62 


1.31 


1353 


1353 


2 


C/1989 Ql Okazaki-Levy-R. 


-12.14± 14.02 


-i-182.62± 14.02 


-h42.90± 


22.24 


-i-80.49± 


27.32 


1.84 


1.30 


458 


452 


2 




-6.47± 12.66 


-i-188.30± 12.66 


-h54.58± 


24.90 


-1-93. 84± 


33.52 


1.86 


1.46 


441 


438 


2 


C/1989 XI Austin 


+22.39± 


2.07 


-365.62± 


2.07 


-h40.84± 


9.89 


-368.96± 


5.42 


1.31 


1.28 


537 


537 


2 




+22.40± 


2.37 


-365.62± 


2.37 


-h41.71± 


12.53 


-373. 15± 


6.77 


1.49 


1.49 


530 


530 


2 


C/1990 Kl Levy 


+17.59± 


1.70 


-857.12± 


1.71 


-1-104.7 1± 


5.72 


-787.37± 


2.88 


1.71 


1.05 


1338 


1323 


3 




+26.02± 


1.65 


-848.67± 


1.65 


-i-111.44± 


6.92 


-790.43± 


3.51 


1.93 


1.21 


1302 


1302 


3 


C/1991 F2 Helin-Lawrence 


+12.53± 


1.79 


-116.77± 


1.79 


-Hl4.24± 


4.03 


-93.14± 


7.31 


0.83 


0.78 


216 


213 


3 




+14.64± 


2.13 


-114.67± 


2.13 


-h24.90± 


4.70 


-89.76± 


12.25 


0.89 


0.86 


212 


212 


3 


C/1993 Al Mueller 


+11.26± 


3.99 


-508.85± 


3.99 


-h62.66± 


2.25 


-408.00± 


2.64 


2.42 


0.98 


1488 


1489 


2 




-19.36± 


3.88 


-539.43± 


3.87 


-h57.95± 


2.74 


-406.83± 


4.16 


2.79 


1.15 


1471 


1471 


2 


C/1993 Ql Mueller 


-i-6.62± 


2.04 


-62.43± 


2.04 


-Hl2.24± 


6.15 


-204.56± 


23.40 


0.97 


0.92 


1047 


1041 


3 




-0.49± 


3.25 


-69.55± 


3.27 


-5.33± 


8.70 


-222.77± 


37.48 


1.21 


1.19 


1039 


1039 


3 


C/1996 El NEAT 


-36.48± 


3.68 


-i-361.78± 


3.67 


-h30.15± 


4.56 


-h368.83± 


5.44 


0.93 


0.60 


495 


492 


2 




-43.35± 


4.05 


-i-354.93± 


4.04 


-Hl8.42± 


6.50 


-h374.29± 


7.89 


1.01 


0.77 


484 


484 


2 


C/1997 BAe Spacewatch 


-i-1.49± 


0.35 


-i-371.77± 


0.35 


-h31.83± 


1.15 


-h402.48± 


1.72 


0.74 


0.67 


1054 


1054 


3 




-1.49± 


0.46 


-i-369.86± 


0.46 


-h38.42± 


1.68 


-h405.01± 


1.62 


1.17 


0.89 


1043 


1043 


3 


C/1997 J2 Meunier-Dupouy 


+38.83± 


0.57 


-2.72± 


0.57 


+44M± 


0.88 


-l-14.72± 


0.91 


0.67 


0.53 


2881 


2863 


2 




+36.86± 


0.57 


-4.69± 


0.57 


-t-45.81± 


1.03 


-l-15.52± 


1.00 


0.74 


0.60 


2824 


2824 


2 


C/1999 Yl LINEAR 


+42.92± 


0.88 


-i-350.42± 


0.88 


-h47.35± 


0.94 


-h345.69± 


1.46 


0.61 


0.48 


1747 


1749 


3 




+19.44± 


0.78 


-i-326.90± 


0.78 


-h49.40± 


1.58 


-i-342.34± 


1.95 


0.92 


0.76 


1737 


1737 


3 


C/2000 SV74 LINEAR 


+50.23± 


0.40 


-85.55± 


0.40 


-h92.31± 


0.85 


-54.88± 


0.60 


1.11 


0.71 


4389 


4349 


3 




+46.66± 


0.39 


-89.12± 


0.39 


-h94.77± 


1.18 


-60.83± 


0.75 


1.23 


0.94 


4356 


4356 


3 


C/2001 Q4 NEAT 


+12.51± 


0.43 


-731.74± 


0.43 


-h60.87± 


0.48 


-696.62± 


0.48 


1.29 


0.63 


5305 


5263 


3 




+05.80± 


0.44 


-738.46± 


0.44 


-h59.59± 


0.70 


-696.74± 


0.64 


1.69 


0.80 


5223 


5223 


3 


C/2002 E2 Snyder-Murak. 


+42.95± 


2.31 


-425.73± 


2.31 


-h89.59± 


12.39 


-441. 69± 


2.79 


0.61 


0.57 


1863 


1863 


2 




+38.11± 


2.79 


-430.57± 


2.79 


-h75.00± 


18.32 


-445. 15± 


3.80 


0.77 


0.75 


1775 


1775 


2 


C/2002 T7 LINEAR 


-13.84± 


0.16 


-653.57± 


0.16 


-h20.72± 


0.39 


-660.08± 


1.05 


0.61 


0.58 


8596 


8768 


3 




-24.42± 


0.24 


-664. 14± 


0.24 


-Hl9.73± 


0.39 


-650.86± 


1.03 


1.21 


0.64 


8643 


8643 


3 


C/2003 K4 LINEAR 


+32.08± 


0.16 


-186.23± 


0.16 


-h30.73± 


0.54 


-126.41± 


0.81 


0.62 


0.54 


7219 


7233 


3 




+40.36± 


0.17 


-177.95± 


0.17 


-h34.73± 


0.63 


-130.71± 


0.77 


0.78 


0.63 


7114 


7114 


3 


C/2004 B 1 LINEAR 


+34.64± 


0.20 


-481.63± 


0.20 


-h35.55± 


0.72 


-460.98± 


0.68 


0.53 


0.43 


3777 


3758 


3 




+34.66± 


0.22 


-481.61± 


0.22 


-h38.01± 


0.81 


-460.07± 


0.81 


0.68 


0.58 


3783 


3783 


3 



specific cases. It appeared, that while qualitative evolution is al- 
most the same, the quantitative results may differ Some examples 
of the significant differences in the output from these two models 
are given in Figs[6]and|7] 

We use here two different frames. The rotating frame is de- 
fined as a heliocentric frame with the i' axis in the radial direction 
pointing toward the Galactic center, the y' axis pointing along the 
Suns velocity, and I' completing a right-handed system. The non- 
rotating frame Oxyz coincides with this rotating one at f = 0, but 



keeps its direction fixed. Basing on iHeisler & Tremaind ( [1983) one 
can write the formula for the tidal force from the disk and the centre 
of the Galaxy in the rotating galactic frame Oxy z as follows 

t = + (A-B) (3A + B)x'x'+ 

-{A-Bfy'y' - [47zGpo - 2{B^ -A^)]z'z' 
where A and B are the Oort galactic constants, G is the gravity 
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Figure 5. Shifts of 1 / a^i i due to the NG effects for 26 of investi- 
gated comets. Three lai'gest uncertainties of l/sLoti.NG belong to comets 
C/1959 Yl and C/1952 Wl and C/1892 Ql (see also Table |2] and sec- 
tion [533t . 




Figure 6. The comparison of the past motion of C/2003 K4 for different 
force-models. The upper (blue) line depicts the swarm of VCs stopped at 
previous perihelion when full Galactic tide (disc + centre) was used. The 
lower (red) line is for simplified (disc only) model. For the comparison the 
starting swarm (very concentrated in this scale) is shown in the right, lower 
corner with black dots. For each swarm the nominal comet orbit is situated 
at the centre of the corresponding dashed-line circle. 



constant, r is the heliocentric distance, p,, is the local disk matter 
density and Mg is the central mass (the Sun and planets), expressed 
in solar masses. Let Gi, G2 and G3 be such that: 

Gi =-{A-B){3A+B) 
Gt,=AkGPo-2{B^ -A^) 

then the equations of motion in rectangular coordinates {Oxyz iner- 
tial frame) are: 




Figure 7. The comparison of the past motion of C/1986 PI for different 
force-models. All symbols have the same meaning as in Fig|6] 



X — G\x cos{Q.ot) + Giy sin(£lof) 
y — G\x! sin(f2of) — G2y cos(f2of) 



if-x 


GMb 






cfly 


GMb 




r3 


dh 


GMb 


dfi 


r3 



3:7 = — r^-z-G3z 



where fl,, is the angular velocity of the Sun about the Galactic cen- 
ter and as a consequence the angular velocity of the rotating frame 
with respect to the nonrotating one and 
.x' = xcos(f2of) + vsin(flof) 
y = — xsin(f2of ) +_vcos(iJo/). 

Substituting this into the equations of motion and remembering 
that G2 = — Gi (all co ntemporary models of Galactic tides adopt 
= A^, see iFouchar d et al. (20(51) as the most recent review of 
the different tide models and their calculation methods) we obtain: 



d^x 


GMb 


1? ~ 


? 


d'-y 


GMb 


dfi 


r3 


dh 


GMb 


dfi 


,3 



' Giz 



which is the set we integrate numerically. We used G\ = 
-7.0702403 X 10"'", and G3 = 4nGp„, where Po, is the local 
Galactic mater density in solar masses per cubic parsec. The model 
limited to the disk tide can be obtained by putting Gi = 0, and 
we used it in parallel, to monitor the separate influence of the 
Galactic centre term. Two examples of such a comparison of dif- 
ferent dynamical models are presented in Figs. [6] and |7] As it con- 
cerns the adopted value of the local Galactic mater density Po, we 
followed most of the contemporary papers in the field, adopting 
Po = 0.100 Mr.,pc^^ . To observe the influence of a possible error 
in this value we repeated the calculations for several comets us- 
ing three significantly different values: Po = 0.050, Po = 0. 100 and 
Po = 0.150 Mp,pc^^. An example of such a comparison is pre- 
sented in Figdl which is significantly magnified to make the differ- 

ence s visible. 

iDybczynskil ( |2006|) have shown, that none of the known stars 
could significantly influence the cometary motion recently (during 
last 10 million years, or so) and in similar time interval in the fu- 
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0.99946 0.99948 0.9995 0.99952 0.99954 

eccentricity 

Figure 8. The comparison of the future motion of C/1996 El with three 
different Po values used. The full Galactic tide model was used here and all 
VCs were stopped at 120 000 AU, going out of the Solar System. The upper 
(blue) swarm was obtained forpo = 0.050 M,.,pc^^ , the middle one (red) 
was obtained for nominal value { = 0. 100 Mcjpc^^) and the lower swarm 
(black one) for Po =0.150 M^pc^^ . The centre of the dashed-line circle 
drawn on top of each of the swarms denotes the position of the nominal 
comet orbit. 

ture. This time interval is comparable with orbital period of an el- 
liptic comet with the semi-major axis of 50 000 AU and thus it is 
close to the timescale of the longest past or fut ure cometary mo- 
tions analyzed in this paper. iDvbczvnskil bOOfih searched for real 
stellar perturbers in all available catalogues, selected several stars 
with maximum possible influence on the Oort cloud comets and 
showed that comparing with Galactic perturbations their influence 
is negligible. He also presented a discussion on the completeness of 
the stellar data (besides positions, both parallax and radial velocity 
are necessary). In short: we cannot be sure that some stars passed 
(or will pass) in the vicinity of the Sun during 20 million years 
time interval, centred at present. But these stars are absent in our 
catalogues e ither because they are too small (e.g. brown d warf, see 
for example: [Research Consortium on Nearby Starsl2009l) or move 
very fast, and now are too far, or both. But only relatively massive 
stars (of order of one solar mass or more) and moving rather slowly 
(less than 30 km s^ ' ) can significantly change cometary orbit, even 
with the aphelion a t 100 000 AU. For more detailed discussion see 
iDvbczviisldl ( |2006|) . 

For the reasons explained above, and taking into account, that 
precise calculation of (even very weak) stellar perturbations would 
be time consuming, we decided to omit them in our model. Per- 
haps in future, especially if some important star will be discov- 
ered (and its parallax and radial velocity will be measured, see 
[Dybczynski & Kwiatkowski 20o3) we will include its gravitational 
influence on the Sun and a comet into our model. 

4.2 Numerical integration 

The original and future orbit calculations produced an output in the 
form of a set of 5001 (nominal orbit plus 5000 VCs) barycentric 
orbits for each comet under consideration. To observe how Galac- 
tic perturbation action varies depending on small initial condition 
changes we separately followed numerically all 5001 VCs, in gen- 
eral, one orbital period to the past and to the futu re. We used the 
well known, widely used and tested routine RA15 ( lEverharlll 19851) 
with the force model described above (with the parameter LL=12 




0.995 0.9955 0.996 0.9965 0.997 0.9975 0.998 



eccentricity 

Figure 9. The past motion of C/2002 T7. Presented here is the distribution 
in q-e plane of all 5001 VCs for that comet, stopped at previous perihelion. 
One can see that the distribution of the previous perihelion distance is very 
close to the Gaussian one. In the center of the dashed-line circle the nominal 
orbit is located. 



to obtain the highest possible accuracy on 64-bit computers). Since 
the Galactic tide is relatively weak, the automatically adjusted in- 
tegration step was generally very large so the speed of the calcu- 
lation allowed us to follow individually such a big number of VCs 
for each comet. 

Since we are interested mainly in the past motion of a comet 
and the apparent source of the observed long-period comets we 
followed their past motion for one orbital period. Having the op- 
portunity to investigate also the future motion of these bodies we 
decided to perform similar calculations for all of them going one 
orbital period to the future. Going further (both back and forth) 
seems useless, since a lot of comets have their previous and/or next 
perihelion distances deep in the Planetary System. Since there are 
no means to calculate precise planetary perturbations for such bod- 
ies say ten million years from now, we cannot follow their motion 
"behind" previous and/or next perihelion point. On the other hand, 
when a comet have large previous/next perihelion distance, their 
subsequent motion do not change our knowledge on the apparent 
source of the observed long-period comets, or their future. 

The situation complicates when past or future orbit of par- 
ticular VC is hyperbolic. We have to stop numerical integration 
at some heliocentric distance, which we call escape limit. To de- 
fine this limit it seems useful to recall contemporary studies on the 
outer border of the Oort cloud. Among rec ent reviews in this field 
we mentio n a chapter i n Com ets II book bv lDones et alj l l2004l) and 
a book bv iFernandej ( |2005|) . Both present consistent conclusions 
that while Oort placed the outer border of the cometary cloud at 
5-15 X 10"* AU current understanding of cometary dynamics makes 
this border much smaller, say 2-8x10'* AU for the aphelion dis- 
tance and l^x 10* AU for semi-major axis. Additionally, the ma- 
jority of the observed Oort spike comets have the semi-major axes 
in the range of 20 000-50 000 AU for pure gravitational orbits but 
this significantly decreases to the value less than 20 000 AU if the 
NG forces are included in the dynamical model. Some discussion 
of this effect is also in Section 5.2 

We decided to be very careful (and a little conservative) so 
we decided to adopt the escape limit to be equal 120000 AU. This 
is the maximum heliocentric distance (or maximum aphelion dis- 
tance) at which we still follow every VCs motion. This was also 
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Figure 10. The past motion of C/ 1989 XI. Presented here is the distribution 
in q-e plane of 4974 returning VCs for that comet (black dots), stopped at 
previous perihelion and 27 escaping (while still elliptic) VCs - dark grey 
(red) dots. One can see that the distribution of the previous perihelion dis- 
tance is significantly different from the Gaussian one. In the center of the 
dashed-line circle the nominal (returning in this case) orbit is located. See 
text for additional explanation. 

applied for very elongated elliptic orbits - if the particular VC 
goes further than 120 000 AU from the Sun we call it escaping and 
stopped the numerical integration. 

This particular choice of the escaping border was unfortunate 
for a handful of comets, where small number of VCs traveled a 
litfle bit further than 120000 AU, but then return to the Sun. We 
examined closely all such cases and if the necessary increase in 
the outer limit seemed acceptable we changed it to obtain homo- 
geneous swarm of returning VCs for that particular comet. Such a 
procedure was applied for example for C/1996 El, where increas- 
ing the outer limit up to 140 000 AU allowed all 5001 VCs to return 
to their previous perihelion among Solar System planets when cal- 
culating its past motion. To obtain the same effect for C/1885 XI 
one should extend the outer limit up to 260 000 AU, what, in the 
light of arguments included in both papers quoted above, seems 
problematic. At such a large hel iocentric distance t he comet mo- 
tion became chaotic, see Fig. 1 in iDvbczvnskil ( 1200 lb . 

Due to the fact, that some VCs swarms consists of both return- 
ing and escaping ones, the rules of stopping the integration were 
sometimes more complicated: the escaping part of the swarm was 
stopped at the 120 000 AU and the rest of VCs was allowed to reach 
previous/next perihelion. There were also a synchronous variant 
of calculation sloping: for returning swarms all VCs were stopped 
when the nominal one reached previous/next perihelion and for es- 
caping or "mixed" swarms, all VCs were stopped when the fastest 
escaping VC reached 120 000 AU. This was done to obtain more 
homogeneous swarm to analyze, where all VCs are stopped at the 
same time. 



5 RESULTS 
5.1 Overview 

Let us start with the context, in which our results should be inter- 
preted. In his paper, lOorll l l 19501) treated original semi-major axes 
as definitive for describing past motion of a comet. In the same 
time he stated, that all comets with 1 /agri < 1 x 10^^ AU^' should 
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Figure 11. The same swarm of VCs for the comet C/1989 XI, but all were 
stopped synchronously when the fastest escaping VC reached 120000 AU. 
Number of returning and escaping VCs is the same of course (4974 and 27, 
respectively). See text for additional explanation. 

be treated as dynamically new, i.e. they are visiting the observ- 
able region for the first time. Such a conclusion came from two 
assumptions. First, the only possible perturbation source are stellar 
perturbations, what means that going to the past we should follow 
original, unperturbed orbit as long as a comet reaches the helio- 
centric distance large enough to be perturbed by passing stars. Sec- 
ond, all observed long-period comets were significantly perturbed 
by stars (if not, their keplerian, original orbits bring them exactly 
to the same perihelion distance at previous passage and no one can 
call them "dynamically new"). 

Today, none of this two assumptions seem to be fully correct. 
As we noted in the Introduction, a comet motion now cannot be 
treated as unperturbed and the dominating (if not the only one) per- 
turbing force is the differential, gravitational action of our Galaxy. 
While Galactic tides perturb the semi-major axis practically on the 
negligible level during one orbital period, the most important part 
of the past and future cometary orbit evolution manifests in the ec- 
centricity changes. 

As a result, while semi-major axis remains almost unchanged, 
the perihelion distance can change significantly, even during one 
orbital period both in past and future motion of course. Taking 
C/1913 Yl as typical example, one can note that the original semi- 
major axis of this comet a= 19022.723 AU and one orbital period 
before it was equal 19022.794 AU. At the same time original per- 
ihelion distance q=1.10AU and one orbital period earlier it was 
equal 2.30 AU. Small changes in the semi-major axis under the 
Galactic perturbations is a well known fact. In the averaged prob- 
lem of the disk tide the semi-m ajor axis is even an integral of 
motion (see iBreiter eTakl (11996') for details). IPretkal d 19981 , 1 19991) 
showed that slow secular changes in the semi-major axis due to the 
action of the Galactic centre can be observed on the much longer 
time intervals than considered here. It means that looking only at 
l/aori of a comet nobody can say what was its previous (or next) 
perihelion distance. 

5.2 Past and future evolution - some statistics 

As it was described in sections |2l [3] and |4l we very carefully de- 
termined orbits for 26 observed long-period comets using sophisti- 
cated observations treatment and NG force model. Than, at osculat- 
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Table 3. The past distributions of tlie NG swarms of VCs in terms of returning [R], escaping [E], including hyperbolic [H] VC numbers; the rule of discrimi- 
nating between [R] and [E] are explained in section l4!2l Aphelion and perihelion distances are described either by a mean value for the normal distributions, 
or three deciles at 10%, 50% (i.e. median), and 90%. In case of mixed swarms the mean values or deciles of Q and q are given for returning part of the VCs 
swarm. For comparison we included the osculating perihelion distance in the second column. Last two columns present the value of 1 / a^n and the percentage 
of VCs, that one can call dynamically new, basing on previous q statistics. 
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a - this pan of swarm includes the nominal orbit, 6-12 VCs have 200 000 AU <Q<320000 AU and 3 VCs - Q between 0.4 - 6.1 x 10^ km c - we used 140 000 AU as the outer border (just 6 VCs 

have Q>120 000 AU). 

ing epoch, we created 5000 VCs for each comet and precisely prop- 
agate all 5001 massless bo dies back and forth in t i me to the limit 
of planetary perturbations. IXodorovic-Juchniewiczl l ll98lh showed 
that the planetary perturbations in the comet's motion should be 
included up to 150-200 AU. Thus, we chose 250 AU as safely dis- 
tant limit where Galactic perturbations are still negligible. Next, we 
followed numerically past and future motion of each comet swarm 
(consisting of nominal orbit and 5000 VCs) under the influence of 
Galactic tides. We stopped the calculation either at past/future per- 
ihelion or when the body in question crosses the assumed "escape" 
border, defined at 120000 AU. 

The results of the investigation of past motion of all 26 comets 
are summarized in Table [3] while the same for future motion is 
shown in Table|4] We present here only results for NG model as we 
are convinced that they gives better description of real cometary or- 
bit and its evolution. For sake of comparison we prepared extended 
version of Tables[3]and|4]consisting of GR and NG results and their 
propagation end-states, but due to the volume of information, this 
is available as the on-line materia l accompanying this paper only 
( iDvbczvnski & Kr61ikowskall2009h . 

Table [3] presents the most important result of this paper. For 
each comet we show, how many VCs return to previous perihe- 
lion (column [R]), how many cross the adopted escape border at 
120000 AU ([E] column) and additionally how many of escap- 
ing VCs have formally hyperbolic orbits at the end of our calcula- 
tion ([H] column, which presents subset of [E] column). Next three 
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Figure 12. The past swarm of VCs for the comet C/1991 F2. Here 2971 re- 
turning VCs (black dots) were stopped at previous perihelion and 2030 es- 
caping VCs (dark-gray or red dots) were stopped when crossed 120000 AU. 
Note that due to the spread distribution of returning VCs, lots of escaping 
VCs are overprinted with each other. 



columns are discussed below. Last column shows the percentage of 
VCs, that one can call dynamically new, basing on previous perihe- 
lion distance statistics, i.e. we calculate the part of dynamically new 
VCs by adding the number of returning VCs with qpiev >15 AU 
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Table 4. The future distributions of the NG swarms of VCs in terms of returning [R], escaping [E], including hyperbolic [H] VC numbers. Aphelion and 
perihelion distances as well as the eccentricity are described either by a mean value for the normal distributions, or three deciles at 10%, 50%, and 90%. In 
case of mixed swarms the mean values or deciles of Q, q and e are given for returning part of the VCs swarm. For comparison we included the osculating 
perihelion distance in the second column. Last column present the value of 1 / ao,i . 
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a - this pait of swarm includes the nominal orbit, h - in this case statistics are for synchronous swarm (see text). 



to the number of escaping VCs. It is remarkable that none of the 
comets from our sample can be called for sure hyperbolic before 
they enter the sphere of visibility. Three of them, namely C/1982 
Ql, C/1952 Wl and 1959 Yl have significant (but always less than 
50 % of all VCs) percentage of hyperbolic VCs and only one, 
C/1952 Wl have formally negative solution for the inverse semi- 
major axis: \/a = (—0.1 ±85.8) x 10^^ AU. However, this value 
was determined with the largest relative uncertainty, what results 
mainly from the smallest number of observations. 

For 15 comets all 5001 VCs are returning and for additional 
5 there is only small (or very small) percentage of escaping VCs, 
constituting typical long tails of a broad distribution. 

For each of 26 comets investigated here we obtain all or sig- 
nificant percentage of returning VCs in their past motion. Thus, 
we were able to present previous aphelion and perihelion distances 
statistics for the returning part of each swarm, see Table[3] In small 
number of cases when the distribution remains Gaussian we present 
its mean value and standard deviation. An illustration of such a case 
is given in Fig. [9] However, in the majority of cases the nonlinear 
action of the Galactic tide have deformed the distributions and they 
cannot be approximated with the Gaussian distribution. In such 
cases we decided to describe these distributions by three deciles, 
namely values separating 10%, 50% (median value) and 90% of all 
VCs for a particular comet. An example of the non-Gaussian dis- 
tribution is given in Fig. [TO] Here 4974 VCs are returning (black 
dots). The strange placement of 27 escaping VCs (dark grey or red 
dots) is caused by the fact, that they were stopped at 120000 AU 
border, while returning VCs were stopped at previous perihelion. 



much later. If we stopped all 5001 VCs synchronously, when the 
fastest escaping VC reached 120 000 AU the distribution of VCs in 
the q-e plane will be quite different, what can be seen in Fig.ll II 

The most complicated cases occur when the swarms of VCs 
divides into two distinct parts: returning and escaping (very often 
with significant number of hyperbolic past orbits). An example of 
such a situation is presented in Figs.ll2land|13l 

As it was already pointed bv lKr61ikows"k3 j200ll.l2004l2006l) . 
taking into account NG forces in the procedure of orbit determina- 
tion shows that LP comets generally have smaller original aphelion 
distances than widely quoted. A typical previous aphelion distance 
varies from 20 000 to 50000 AU, what makes the average semi- 
major axis to be less than 20 000 AU. This is a significantly smaller 
value that typically thought for so called Oort cloud comets. 

However, the most important information about the prove- 
nance of long-period comets is their previous perihelion distance. 
It is really striking that 50% of comets in our sample definitely vis- 
ited the Planetary System at their previous perihelion, sometimes 
as deep as I AU (or less) from the Sun. No one should again call 
them dynamically new comets! One of the most remarkable ex- 
ample is comet C/1990 Kl, commonly considered as dynamically 
new and discussed in detail in Section |5.3.6| Additionally, for the 
remaining half of our sample, only six comets, namely C/1946 UI, 
C/1956 Rl, C/1991 F2, C/1993 Ql, C/1997 BAf, and C/2002 T7 
seem to be dynamically new with high degree of certainty. Re- 
maining 7 comets have mixed past swarms of VCs, with the per- 
centage of the previous perihelion distance greater then 15 AU (the 
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Figure 13. The same swarm of VCs for the comet C/1991 F2 as in Fig[T2l 
but all were stopped synchronously when crossed 120000 AU. 
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Figure 14. The swarm of VCs for the future motion of C/1952 stopped 
synchronously when the fastest escaping VC crossed 120000 AU. 

widely adopted threshold value for significant planetary perturba- 
tions) varying from 25% to 63%. 

Performing a very simple (while really impressive) statistics, 
one can sum all VCs for all 26 comets and calculate the percentage 
of VCs with previous perihelion distance larger than 15 AU. The 
result is 32.9%. Statistically this is the percentage of dynamically 
new comets among observed long-period comets, basing on our 
sample of 26 comets. 

This value i s sig nificantly smaller than obtained by 
iDvbczviiskil ( 1200 ll : l2006h . who claimed that about 50% of long- 
period comets with semi-major axes larger than 10000 AU could 
be treated as dynamically new. This discrepancy shows the im- 
portance of NG forces applied at the stage of orbit determination 
jPvbczviiski 2001, 2006 used pure gravitational orbits). It is re- 
ally impressive if we keep in mind that gravitational original orbit 
of eight comets from our sample are hyperbolic and only four of 
them might be called dynamically new when NG original orbit is 
properly accounted for. 

As it was mentioned earlier, we also studied the future motion 
of all 26 comets in our sample, in an analogous manner. The re- 
sults of this part of investigation are summarized in Table |4l The 
description is much simpler here, because we have only one comet, 
C/1952 Wl, with significantly mixed swarm of VCs, consisting 




4.3 -i r 

0.99992 0.999925 0.99993 0.999935 0.99994 0.999945 0.99995 
eccentricity 

Figure 15. The swarm of VCs for the future motion of C/1997 J2 stopped 
synchronously when the fastest escaping VC crossed 1 20 000 AU. 

of 1509 returning VCs and 3492 escaping (with 3229 hyperbolas 
among them, including the nominal orbit). A synchronous plot of 
all 5001 VCs for this comet is depicted in Fig[T4] It seems that this 
comet will be lost in the future. The similar fate can be attributed 
for next 16 comets having 100% of escaping, hyperbolic VCs and 
for C/1885 XI (over 99% escaping hyperbofic VCs). 

Only 6 comets, C/1913 Yl, C/1946 Ul, C/1986 PI, 
C/ 1996 El, C/1997 BAe and C/1999 Yl have 100% of returning 
VCs and all of them will have the future perihelion distance deep 
among Solar System planets. The future motion of these comets 
seems to be unpredictable due to possible planetary perturbations. 

Two special cases of future motion: C/1989 Ql and C/1997 J2 
are discussed in the next section. 

Summarizing, from our sample of 26 comets 17 should be 
treated as lost in the future, one (C/1997 J2) may remain bound 
to the Solar System on large elliptical orbit and 8 will return in the 
future, pretending they are coming from the so called inner Oort 
Cloud, and having perihelion distances well inside the sphere of 
visibility. 



5.3 Individual cases 

5.3.1 C/1885 XI Fabry 

This is the oldest comet investigated in this paper. Its past mo- 
tion is almost clear. When dealing with different 5000 VC swarms 
of this comet, sometimes we obtain just few (or a dozen) of 
VCs with hyperbolic past orbit. In every case, there is also some, 
very small, percentage of VCs, which have their aphelia slightly 
above the outer limit of the heliocentric distance, adopted in this 
paper as the limit for returning comet (120000 AU). The rea- 
son for this is the fact, that the A2 non-gravitational parameter 
is hard to determine and have relatively large uncertainty (A2 = 
(-0.2879 ±0.1606) x 10"** AUday^^ ). Together with relatively 
large mean error (3.58 arc-sec, anyway, very good for such an old 
comet!) it causes the VC swarms to be more spread than in many 
other cases. Assuming, that the A2 value is still acceptable, we per- 
formed two additional investigations for this comet, using two non- 
standard VC swarms: first produced with fixed A2 value and the 
second, with fixed both A; and A2. The result of this additional cal- 
culation is very impressive, what is clearly visible in Fig[T6] The 
swarm obtained with fixed A2 value appeared much more compact. 
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Figure 16. The swarm of VCs for the pas motion of C/1885 XI stopped at 
previous perihelion. VCs were generated with the non-gravitational param- 
eter A2kept constant. 



all VCs are returning in their past motion even for standard limit of 
120000 AU. Additional fixing of value A\ in the second case did 
not change the situation significantly (it has much smaller error of 
determination). For example, the previous perihelion distance for 
standard swarm (17 VCs escaping) varied from near zero to almost 
1000 AU (however 90% of VCs have q < 4.3AU). For the swarm 
obtained with fixed A2 value all VCs are returning and have the 
previous perihelion distance less than 120 AU, with only 3 of them 
greater than 15 AU and 90% smaller than 1.72 AU! 

Taking all this into account, we can state with high degree of 
certainty that this comet was not a dynamically new one, having the 
l/a^-j = (61.0± 17.2)x lO^^AU^' and the most probable previous 
perihelion distance of about 1-2 AU, i.e. deep among planets. It 
may be worth to mention here, that this is one of the famous 19 
long-period comets which Oort used to support his hypothesis in 
1950. 

The future motion of this comet is much simpler to describe. 
The vast majority (99.9%) of VCs in our standard model are es- 
caping in their future motion (almost all of these - 99.3% of 
all VCs - are hyperbolic). The single, smallest returning orbit 
has qfut = 14AU and Qfut =60000AU. The l/fl/,„ = (-102.0± 
41.1)x lO^^AU^' here, so this comet is almost certainly lost from 
the Solar System. 



5.3.2 C/1913 Yl Delavan 

Previous perihelion passage of this comet is determined with the 
highest accuracy of a group of comets discovered before 1970. 
Our calculations show that C/1913 Yl was previously inside in- 
ner Solar System and lived it in an orbit of qprev <3 AU and 
Qprev <43 000 AU (Table O This is the second comet in our sam- 
ple which was originally considered by Oort as dynamically new 
object. It seems worth to mention here, that such a small previous 
perihelion distance was obtained as a result of non-gravitational 
forces incorporated in this paper for the original and future orbit 
determination. Pure gravitational model gives the previous perihe- 
lion value of the order of 150 AU for this comet. 

The future motion of C/1913 Yl is also determinable with 
very high accuracy. In contrary to C/1885 XI here all VCs are re- 
turning with the distribution of the next perihelion and aphelion dis- 
tances only slightly departing from the Gaussian one. However, by 
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Figure 17. Past and future motion of all 5001 VCs for C/1913 Yl in terms 
of their heliocentric distance. Solid, central black line depicts the nominal 
orbit and black histograms show the density of VCs previous/next perihe- 
lion passages in time. The time axis is centered (t=0) at the observed peri- 
helion passage of this comet in 1914. 



fitting to Gaussian one can derive the values of 0.97 AU ±0.02 AU 
and (23.1±l.l)x lO-' AU for next perihelion and aphelion distance, 
respectively (compare with Tablel?^. The heliocentric motion of all 
VCs for one past and one future orbital periods is shown in Fig|17l 
If we demand, that the Oort cloud comet should move outside 
the sphere of planetary perturbations (before it is observed for the 
first time) C/1913 Yl do not fulfill this condition either in the past 
or in future. 



5.3.3 Comets C/1892 Ql. C/1952 Wl and C/1959 Yl 

These three old comets have the past swarms of VCs widely spread. 
Apart from a large number of escaping VCs, the significant part of 
them have hyperbolic orbits. While in all cases we can call them 
'dynamically new' with rather large degree of certainty, the exact 
provenience of these comets is hard to be determined, we cannot 
exclude even their interstellar origin. The probability of such a case 
was extensively discussed bylHughes ( 1991). 
C/1952 Wl is the only comet in our sample which have for- 
mal value of l/flori < (see Table O. The future motion of 
C/1892 Ql and C/1959 Yl show evident hyperbolic ejection while 
for C/1952 Wl 64% of VCs are escaping on hyperbolic orbits but 
on the other hand some 30% of VCs will return, sometimes with 
very small perihelion distances (see Tablel?^. 



5.3.4 C/1 986 PI Wilson 

This is very difficult comet to investigate because of its split- 
ting into two fragments. According to MWC08, the main compo- 
nent A of this comet exhibits NG acceleration in the time interval 
1986 August 05 - 1989 April 1 1 . Therefore, we also determined the 
NG effects from the whole observational arc (see TablefTJ and next 
starting from these nominal orbits (GR and NG) the past and fu- 
ture orbital evolution have been analyzed. It turns out that the frag- 
ment B was a small piece of parent body and had ceased all activity 
or had disintegrated after April 1989. To estimate the possible un- 
certainty connected with this splitting we tried to find the nominal 
orbits from two separate sets of data: before and after splitting, but 
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Figure 18. Future swarni forC/1986 PI. Note how small differences appeal' 
between VCs during the next perihelion passage roughly 55 000 year from 
now. The centre of the dashed-line circle marks the nominal orbit. See text 
for additional comments. 

the exact date of this break-up is not clear. During comet's obser- 
vation on 1988 Fe bruary 13 wa s disco vered that nucleus was frag- 
mented. However. iMeech et al J i ll 99^ argued that the most prefer- 
able date for splitting was mid-October 1987 because of the likely 
connection with a brightness of outburst in that time. Thus, we de- 
cided to determine the past orbit of C/1986 PI from the time in- 
terval 1986 August 05 - 1987 July 04 (next astrometric observation 
was taken three months later on October 2) and the future orbit 
from the time interval 1987 October 26 - 1989 April 11 (114 ob- 
servations). We derived values of (53. 14 ±4. 14) x lO"*" AU"' and 
(30.12± 1.84) X lO"*" AU^i for l/aori.NG and l/aoiicR, respec- 
tively, where NG solution was determined with rms= 1 "30 and pure 
gravitational solution - with rms=l"35. In this approach we de- 
rived 37.6l;^ { X 10^ AU and 1.90lg-^| AU for previous aphelion 
and perihelion passage, respectively, where by adding the lower 
and upper errors to nominal value the deciles of 10% and 90% 
can be derived and directly compared with respective values in Ta- 
ble [3] By comparing both results of previous perihelion swarms 
we stated that the comet C/1986 PI were deeply inside inner solar 
system with perihelion not larger than 4 AU and aphelion less than 
5 X lO'* AU. 

For the data taken after the splitting the NG effects were inde- 
terminable. The VCs swarm of future motion is the most compact 
in our sample, what can be observed in Fig[T8l We derived the value 
of (749.48 ±9.83) x 10"^ AU"' for l/afut,GR By comparing this 
value with result given in Table |4] we conclude that this comet will 
start subsequent perihelion passage with small q (approximately 
1.2 AU) and with aphelion distance Q~3 000AU, very similar to 
future aphelion distance of Comet C/1996 B2 Hyakutake. 

J..?. J C/1989 Ql Okazaki-Levy-Rudenko and C/1989 XI Austin 

Both comets were observed in short time intervals (four and six 
months, respectively) thus the past and future perihelion swarms of 
VCs are relatively disperse what gives the formal chance of only 
about 25% that these comets are dynamically new. 
Comet C/1989 Ql have a great majority (99%) of returning 
VCs with previous perihelion distance of the order of 1 AU but 
there is a long tail in its distribution, consisting of 51 escaping 
VCs, including 6 hyperbolas. Thus, the resulting ^-deciles are 



0.41 AU-0.55 AU-33.4AU (at 10%, 50%, and 90%) and taking 
limit ^previous =15 AU for dynamically new comets the probabil- 
ity of 24.4% was formally obtained that this comet is dynamically 
new. 

Similarly, the q deciles of 1.7 AU-5.9 AU^9.4 AU (at 10%, 50%, 
and 90%) give the chance of 25.9% that comet C/1989 XI is dy- 
namically new. 



5.3.6 C/mOKlLevy 

This comet may serve as an excellent illustration, how the incorpo- 
ration of NG forces into the process of orbit determination shrinks 
the semi-major axis and significantly shorten the time interval of 
the past motion. When we use pure g ravitational solution for this 
come t (see the on-line extended table. iDvbczviiski & Krolikowskal 
|2009|) only 3644 VCs are classified as returning while 1357 as es- 
caping but elliptical. For the slowest returning VC it takes almost 
15x 10^ years to reach the previous perihelion. If we want all VCs 
to become returning, the 'escape limit' should be increased up to 
190000 AU what seems unacceptable. But it is worth to mention 
that in that case the slowest VC takes over 27 million years to pass 
the previous perihelion. 

i^From Table[3]one can read, that when we use NG model, all 
VCs are returning with the standard 'escape border'. Moreover, the 
slowest VC pass the previous perihelion point just 1.25 x 10* years 
ago! 



5.3. 7 C/1991 F2 Helin-Lawrence and C/1993 Ql Mueller 

These comets seem to be classical jumpers from the Oort Cloud 
with largest Aq = qpievious ^ qori among considered objects (larger 
than 300 AU). 

Since the astrometric set of data for C/1991 F2 consists of 
just 1 14 observations this comet represents a limit case for weight- 
ing. We decided to show the past and future motion on the basis 
of weighted observations (Tables [3] and |4l(. It should be noted how- 
ever, that in the non-weighted data case the past swarm of VCs are 
almost all returning (4968 of 5001 VCs) with q range: 22.3 AU- 
71.9AU-340.8 AU (three deciles at 10%, 50%, and 90%) what 
gives 98.8% of dynamically new VCs. Thus, the main conclusion 
that this comet is dynamically new remains unchanged, indepen- 
dently of data treating. 

The past swarm of comet C/1993 Ql is dominated by escaping 
VCs with maximum of aphelion distance ~150000 AU (compare 
with returning part of swarm given in Table[3](. Our results suggest 
that this comet at previous perihelion was extremely far from the 
Sun. 

Both comets will be probably lost in the future - they leave 
the solar system on hyperbolic orbits. 



5.3.8 C/1996 El NEAT 

Because only 6 VCs of this comet crossed the adopted escape limit 
in the past motion we decided to increase this limit up to 140000 
AU when performing backward numerical integration. Then, all 
5001 past VCs were returning and 74% of them had q<15 AU. 
After the observed passage through the Solar System its nominal 
semi-major axis is drastically reduced from 54000 AU down to 
2600 AU and it will pass the perihelion 1.3 AU from the Sun. 
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5.3.9 Comets with large observed perihelion distance 

In our sample there are four comets with relatively large osculating 
perihelion distances: C/1997 BAe Spacewatch, C/1997 J2 Meunier- 
Dupouy, C/1999 Yl LINEAR and C/2000 SV74 LINEAR. First 
three might be grouped together to show different dynamical his- 
tories of these comets in spite of their very similar original perihe- 
lion distances (from 3.05 AU to 3.44 AU) and rather similar original 
semi-major axes (within the range from 21 000 AU to 31 500 AU, 
Tab|2](. One can see in Tables [3}@] that the previous and future per- 
ihelion passages are well determined for all of them. Thus, we can 
state that only the first one (C/1997 BAg Spacewatch) having pre- 
vious perihelion distance around 20.0AU±3.3AU seems be dy- 
namically new, whereas remaining two comets were deep in the 
inner Solar System during the past perihelion passage (Table [Sj. 
This may serve as an additional argument, that even having all pa- 
rameters of the original orbit one cannot judge what was the previ- 
ous aphelion or perihelion distance of a particular comet. While the 
strength of the Galactic perturbations depends mainly on the semi- 
major axis value, it can be significantly different, depending on the 
orientation of the orbit with respect to the Galactic equator. There 
exist another point, from which C/1997 J2 Meunier-Dupouy is also 
an interesting comet. Its previous perihelion distance (2.80 AU) 
was slightly smaller than the original one (3.05 AU), what means 
that this comet was discovered during the "growing" phase of its 
perihelion distance evolution under the influence of the Galactic 
tides (two more comets in our sample, C/1978 HI and C/1990 Kl 
have also smaller previous perihelion distance than the observed 
one). The statistical analysis of the q evolutional phase of long- 
pe riod comets at their discove ry epoch was performed for example 
by iMatese & Lissaue3 ( l2004h . They concluded that most comets 
should be observed during the "decreasing" phase but noted sev- 
eral opposite situations. It is worth to mention however, that their 
analysis was based on pure gravitational cometary original orbits. 
Future history of C/1997 J2 is more uncertain. It can be seen from 
Table |4] that next VCs aphelion distances are slightly larger than 
our assumed standard outer border: 97% of VCs are classified as 
escaping, but formally all of them are elliptic, see Fig[T5] By shift- 
ing the escape border to 170 000 AU we can obtain all VCs return- 
ing and approximately Gaussian distribution of aphelion distance 
(2 = 136.3±8.4 thousand AU) and non-Gaussian for perihelion 
distance with values of 5.6AU, 7.9 AU, 15.9 AU for P', 5* and 
9* decile of q, respectively. But so large aphelion distance means 
that this comet will probably leave the Solar System. As it can be 
noted from Table |4l two remaining comets, C/1997 BAg Space- 
watch and C/1999 Yl LINEAR, will have a very similar future 
evolution. Both of them in the subsequent perihelion passage will 
enter the inner Solar System in orbits with perihelion distances be- 
tween 3 AU and 4 AU and aphelion of ~5 000 AU-6 000 AU. 

Fourth comet with large osculating perihelion distance, 
C/2000 SV74 LINEAR, have significantly shorter original semi- 
major axis of about 10 800 AU and its previous perihelion passage 
was below 4AU. However, it seems that C/2000 SV74 LINEAR 
will leave the solar system on a hyperbolic orbit. 



5.3.10 C/2001 Q4 NEAT and C/2002 T7 LINEAR 

These comets were at great distance from the Sun when discovered 
(comet C/2001 T7: over II AU, C/2002 T7: ~7AU) and were al- 
most stellar in appearance. At maximum brightness both became 
easy binocular objects (4 mag. and 2.5-3 mag, respectively) and as 
consequence were observed long after perihelion passages. Such a 



long observational arcs allowed to determine their orbits with high- 
est precision and we obtained very accurate results with Gaussian 
distribution of VCs parameters both for past and future motion. 

Both comets have past swarms of VCs fully returning and 
compact, however they differ dramatically in previous perihelion 
distance. Comet C/2001 Q4 was at previous perihelion close to 
the Sun (1 .87 AU±0.06 AU), whereas CI 2002 T7 passed the previ- 
ous perihelion at a distance of 142 AU±19 AU from the Sun. Thus, 
from this pair, only C/2002 T7 seems to be dynamically new. 

While for these comets different form of the g{r) functions 
might be more appropriate, we are convincec^] that the main con- 
clusions remains the same. 

Both comets are leaving the Solar System on hyperbolic or- 
bits. 



6 CONCLUSIONS 

The aim of our research is to refine widely quoted opinions about 
the apparent source of the long period comets and a popular habit 
to call some of them dynamically new by looking on the original 
inverse semi-major axis value only. 

In this paper, the first one of a planned series, we touch 
the most important group of comets, 26 long-period comets with 
1/aoi-i < lO^'^AU^i for which we can obtain very precise or- 
bits including NG forces in the dynamical model. The impor- 
tance of this sample comes from the fact that for almost all of 
these comets we obtained significantly smaller semi-major axes 
due to NG force incorporation in the process of orbit determina- 
tion, but still all of them constitute the NG Gort spike, i.e have 
1/aori.NG < lO-^AU-l. 

Having such a precise orbits we applied ingenious Sitarski's 

( 1 19981) method of creating swarms of VCs, all compatible with ob- 
servations. This allowed us first to obtain original and future or- 
bits parameters with their uncertainties and then to propagate all 
individual VCs one orbital period to the past and to the future, by 
means of strict numerical integration of the equation of motion with 
the Galactic tide action included. 

Note, that our neglectin g of stellar perturb ations is based on 
well documented arguments jPybczy hski 2006), but the incorpo- 
ration of NG forces into the orbit determination makes these ar- 
guments even str onger. All but one o f comets in our sample (and 
many others, see lKr61ikowskall2006l) appear to have significantly 
smaller semi-major axes what implies that their orbital period be- 
comes drastically shorter. This makes a strong stellar perturbation 
by an undiscovered massive and slow moving star extremely im- 
probable. 

Our results show, that two widely quoted opinion about the 
long-period comets are not necessary the "canonical truth". First, 
all 26 comets analyzed in this paper have gravitational inverse of 
the semi-major axis very small (l/aori,GR < 6 x 10^^ AU^') what 
makes almost all authors to call them dynamically new. But from 
Table[3]one can clearly recognize, that no more than 5 or 6 of them 
deserve for such a name. The rest visited our Planetary System dur- 
ing the previous perihelion passage, possibly experiencing strong 



^ We have tested various shapes of others g(r)-like functions (Eq[2) for 
both comets and differences in comet's osculating orbital elements among 
all considered NG models were about one order of magnitude smaller than 
obtained between gravitational and NG models 
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planetary perturbations and nobody knows what process or mech- 
anism direct them at that time (or later) to the observable orbit at 
which we discovered them. 

Second, we showed, that about 30% of our sample of long- 
period comets will remain observable in the next perihelion pas- 
sage. Assuming that this is not a strange statistical fluctuation and 
keeping in mind that the only recognized force which makes their 
orbits to evolve is the tidal action of the Galactic matter, the widely 
quoted necessary condition for a comet to be made observable by 
Galactic perturbation is too strong. It is typical to demand that such 
a comet should decrease its perihelion distance from above 15 AU 
down to the observational value in one orbital period, otherwise 
they cannot cross the so called Jupiter- Saturn barrier, what now 
seems to be a "not necessary" condition: 

We have shown that a significant percentage of long-period 
comets can visit the zone of visibility during at least two or three 
consecutive perihelion passages. 

In this context we would like to suggest some correction to 
the widel y quoted and ver y useful comet classification scheme pro- 
posed by iLevisonI l ll996l) . In his taxonomy the nearly isotropic 
comets (NICs) are divided into two subpopulations: new NICs and 
returning NICs on a basis of their original semi-major axis. This 
division should probably be replaced by a new one, based on the 
past dynamical history of a particular comet, especially if we want 
to keep naming NICs 'new' or 'returning' from the dynamical point 
of view. 

There exist also a difficult problem of linking the 'dynamical 
age' of long-period comets with their physical properties, includin: 
the chemical com position jCrovisie 3 l2007l : ICrovisier et al.|[200^ 



iBiver et al We tried to look for such a linkage but it seems 

that physical and dynamical characteristics might be quite different. 
For example the spectra of C/2001 Q4 NEAT and C/2002 T7 LIN- 
EAR were intensively investigated as they were bright and well ob- 
served comets/nieobservations of large molecules in both comets 
show ( iRemiian et alj|2006h that they are chemically very similar 
while our 'dyna mical ages' for these comets are completely dif- 
ferent. However, iRemiian et al.l j2006l) also stated that when com- 
paring the molecular production ratios with respect to water it ap- 
pears that C/2002 T7 is more similar to Hale-Bopp physical class of 
comets (comets rich in HCN: HCN/H20> 0.2%), while C/2001 Q4 
is more similar to Hyakutake class (comets with HCN/H2O abun- 
dance ratio of ~0.I%) where two more investigated here comets 
(C/I989 XI and C/I990 KI) are also included. On the other hand, 
C/1995 01 Hale-Bopp and C/1996 B2 Hyakutake seems to be dy- 
namically more similar than comets Hale-Bopp and C/2002 T7. 

^From the point of view of the apparent source of the long 
period comets, the main result of our paper is the list of previous 
aphelion and perihelion distances together with numbers of escap- 
ing VCs, presented in Table [3] But, comparing the past and future 
(Table |4ll dynamics of 26 comets bring us to an interesting addi- 
tional conclusion. 

Treating these results as a probe of the past and future flux 
of observable long-period comets one can conclude that there ex- 
ist surprising discrepancy between past and future cometary orbit 
evolution. Some 33% of comets arriving at the outer border of the 
Planetary System (250 AU from the Sun) are dynamically new, thus 
66% have visited our Planetary System one orbital period before. 
On the other hand, 65% of comets leaving the zone of visibility will 
be lost in interstellar space and only 35% will return. But, in the 
"next turn", they should constitute 66% of observable long-period 
comets, what might suggest that the flux of long-period comets 
should constantly decrease... 



And the last remark, of more histori cal than scientific value. 
As far as we are able to recognize now l lDvbczvnskj||2009l) only 
two comets from our sample studied in this paper were used by 
lOorj tl950j in construction of his Table 1, namely C/1885 XI 
Fabry and C/I913 Yl Delavan. According to our results both these 
comets have visited the Planetary System during their previous per- 
ihelion and passed rather close to the Sun (1-2 AU) approximately 
2.5 X 10^ years ago. 

As it concerns our future plans in this research, we have cal- 
culated over 100 precise original and future orbits, waiting for the 
detailed analysis of their past and future motion. In the next step 
we probably concentrate on "statistics improving" by significant 
increasing of number of analyzed comets, with a special attention 
on large perihelion distance comets, where even without detectable 
NG effects, orbits are very precise. 



7 APPENDIXES 

7.1 Chauvenet's and Bessel criterions 

Both criteria differ in the upper limit of the accepted residu- 
als, e.g. observed minus computed values of right ascension, 
Aa ■ c os S, and declinati on, A8. According to the Chauvenet's cri- 
terion Ich auvenetllT908h from the set of A' residuals, we should 
discard all values of ^ for which 



\^\>a-Kij2{N) 

where tr is a dispersion of ^ : 



\ 



Z^k2 I /N 



and /2 {N) is the unknown upper limit of the integral of the prob- 
ability distribution, ^{t,): 



(i>{x)6x =1 , 

where x = t, /a. 

According to this criterion the data point is rejected if the 
probability of obtaining the particular deviation of residuals from 
the mean value is less than I/(2A'). To determine this probability 
the normal distribution of <^ is assumed. 

The Bessel criterion (more restrictive than the Chauvenet's crite- 
rion) rejects from the set of A' residuals all the values of t, for which 

\£,\>a-Kx{N), 

where K\ [N) is defined by: 



1 

N' 
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